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Atty.DktNo.:301-0005US 
NMR SPECTROSCOPY USING SPIN STATES WITH LONG LIFETIMES 

CROSS-REFERENCE TO RELATED APPLICATIONS 
[000 1 ] This application is a National Stage of PCT/GB2004/002966, entitled "NMR 
Spectroscopy Using Spin States with Long Lifetimes" and having an international filing date 
of 8 July 2004, to which priority is claimed and which claims priority to GB03 17600.5 (28 
July 2003) and GB0329692.8 (22 December 2003). 

FIELD OF INVENTION 
[0002] The present invention relates to the field of Nuclear magnetic resonance 
(NMR). More specifically, the present invention relates to a method of Nuclear magnetic 
resonance (NMR) spectroscopy which utilises non-equilibrium spin states with long 
lifetimes. 

BACKGROUND OF INVENTION 
[0003] Nuclear magnetic resonance (NMR) is a powerful and well established tool 
for studying samples and sample interactions. It can also be used for imaging, for example, 
Magnetic Resonance Imaging (MRI). 

[0004] In NMR, the spin and magnetism of atomic nuclei are exploited to provide 
information about the chemical composition, spatial distribution, or molecular motion of the 
molecules or atoms. 

[0005] One of the limitations of NMR is the low intrinsic signal strength. Attempts to 
address this problem have included nuclear spin hyperpolarization experiments, for example 
using parahydrogen (C. R. Bowers and D. P. Weitekamp, Phys. Rev. Lett. 57, 2645-2648 
(1986); C. R. Bowers and D. P. Weitekamp, J. Am. Chem. Soc. 109, 5541-5542 (1987); M. 
G. Pravica and D. P. Weitekamp, Chem. Phys. Lett. 145, 255-258 (1988)). 
[0006] One of the principal factors contributing to the applicability of NMR based 
techniques is the relatively long lifetime of the nuclear spin order. This long lifetime allows 
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NMR to be used to follow processes such as diffusion, flow, slow molecular motion and 
chemical reactions. The relaxation of the nuclear spins back to thermal equilibrium, and 
hence the lifetime of the nuclear spin memory is characterised by a time constant, Ti known 
as the longitudinal relaxation time constant or as the spin lattice relaxation time constant. 
The time constant Tj is governed by molecular motion, and is in general different from the 
phase memory time T2, which depends upon molecular motion and also external 
perturbations such as magnetic field gradients. In general, Ti is always longer than O.5T2, 
although in common situations, Tj is much greater than T 2 , often by an order of magnitude. 
[0007] Generally, for many samples, the time constant T\ is of the order of seconds. 
Thus, the time for following a process of the type mentioned above is very limited. 
[0008] Attempts have been made to extend the phase memory time in the context of 
quantum computation, using the concept of Decoherence-Free Subspaces (L. Viola, E. M. 
Fortunato, M. A. Pravia, E. Knill, R. Laflamme and D. G. Cory, Science 293, 2059-2063 
(2001); E. M. Fortunato, L. Viola, M. A. Pravia, E. Knill, R. Laflamme, T. F. Havel and D. 
G. Cory, quant-ph /0210057 (2002); D. A. Lidar and K. B. Whaley, quant-ph /0301032 
(2003); J. E. Ollerenshaw, D. A. Lidar and L. E. Kay, quant-ph /0302175 (2003)). However, 
although the compensation of artificial noise contributions to the phase memory time T 2 has 
been demonstrated, the extension of the spin memory time beyond Ti has not been 
demonstrated. 

SUMMARY OF INVENTION 
[0009] The present inventors have found that it is possible to generate nuclear spin 
states which remain far from equilibrium for a time much longer than T\. These spin states 
may be advantageously used to study processes with a longer lifetime than that of Tj. These 
non-equilibrium spin states may be studied (and in some cases, created) by switching the 
symmetry of the nuclear spin Hamiltonian operator. 

[00 1 0] Thus, in a first aspect, the present invention provides an NMR method, the 
method comprising: providing a sample where the nuclear spin Hamiltonian operator of the 
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component molecules of the sample possess one or more symmetry operations; creating a 
quasi equilibrium nuclear spin ensemble state in a sample, said quasi equilibrium nuclear 
spin ensemble state comprising at least two manifolds of spin states which transform 
differently under said symmetry operations of said Hamiltonian and said manifolds having 
different mean nuclear spin populations; allowing said quasi equilibrium nuclear spin 
ensemble state to remain for a time of equal to or substantially greater than 3Ti, where Ti is 
the spin lattice relaxation time; breaking at least one of the said symmetry operations of said 
Hamiltonian; applying a sequence of magnetic fields to generate a nuclear magnetic 
resonance signal from said sample; and detecting said nuclear magnetic resonance signal. 
[001 1 ] The term symmetry is used here to describe an exchange operation of one or 
more of the nuclear spins that leaves the nuclear spin Hamiltonian invariant. In practice, the 
imposition of a new symmetry operation does not need to be exact. This is discussed in 
more detail later. It will be appreciated by those skilled in the art that the said Hamiltonian is 
the effective Hamiltonian of the nuclear spins, which governs the evolution of the nuclear 
magnetic moments over the relevant part of the applied experimental sequence. 
[001 2] Manifold is the term used in the art for groups of nuclear spin states that 
transform in a common way under the symmetry operations of the nuclear spin Hamiltonian. 
The manifolds may reach internal equilibrium independently and may have different mean 
nuclear spin populations. A state of internal equilibrium for the individual manifolds but 
different mean nuclear spin populations for different manifolds is, by definition, not a 
thermal equilibrium state for the entire ensemble of nuclear spins. Nevertheless, it may be 
very long lived, since there are no relaxation processes internal to the molecules which can 
lead to its destruction, hence it is referred to as a quasi equilibrium nuclear spin ensemble 
state. 

[001 3] For the avoidance of doubt and as used herein, the term quasi equilibrium 
nuclear spin ensemble state will be used to mean a state of the spin ensemble in which 
manifolds of different symmetry have reached internal equilibrium but in which the different 
manifolds are not equilibrated with each other. 
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[00 1 4] An example of a quasi-equilibrium spin state is found in an ensemble of 
molecules each containing two spins- 1/2 of the same isotopic type, in a magnetic field which 
is sufficiently low that the chemical shift differences between the nuclei may be ignored. In 
this case, an exchange of the two nuclear spins is as a symmetry operation of the nuclear spin 
Hamiltonian operator. The nuclear spin states may be classified as symmetric with respect to 
this operation (the three triplet states) or as antisymmetric with respect to this operation (the 
singlet state). A state of the ensemble of molecules, in which the populations of the triplet 
states have equilibrated with each other, but in which the singlet population is different from 
the mean triplet population, qualifies as a "quasi-equilibrium nuclear spin state". 
[001 5] Another example is found in a sample composed of molecules of the form A- 
CH2-CH2-X, where A and X are arbitrary groups, supposed to have no significant nuclear 
spin magnetism. In low magnetic field, an exchange of the spins in the two CH 2 groups with 
respect to each other is a symmetry operation of the nuclear spin Hamiltonian operator. The 
nuclear spin states may be classified as symmetric, or antisymmetric, with respect to this 
spin exchange operation. A state of the ensemble of molecules, in which the populations of 
the symmetric states have equilibrated with each other, and in which the populations of the 
antisymmetric states have equilibrated with each other, but in which the mean symmetric 
state population is different from the mean antisymmetric state population, qualifies as a 
"quasi-equilibrium nuclear spin ensemble state". 

[001 6] A non-equilibrium singlet population is a special case of a quasi-equilibrium 
spin state, which applies when the molecules only contain two nuclei with significant 
magnetic properties. Quasi-equilibrium nuclear spin ensemble states apply to molecules 
with two or more coupled nuclear spins. The present invention may be used with a sample 
which can form a non-equilibrium singlet population or with a sample which can form a 
quasi-equilibrium spin ensemble state other than a non-equilibrium singlet population. 
[001 7] In a second aspect, the present invention provides an NMR method, the 
method comprising: providing a sample where the nuclear spin Hamiltonian operator of the 
component molecules of the sample possess one or more symmetry operations; creating a 
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quasi equilibrium nuclear spin ensemble state in a sample, said quasi equilibrium nuclear 
spin ensemble state comprising at least two manifolds of spin states which transform 
differently under said symmetry operations of said Hamiltonian and said manifolds having 
different mean nuclear spin populations, said quasi equilibrium nuclear spin ensemble state 
not having a singlet state population; breaking the symmetry operation of said Hamiltonian; 
applying a sequence of magnetic fields to generate a nuclear magnetic resonance signal from 
said sample; and detecting said nuclear magnetic resonance signal. 
[001 8] In the method of the second aspect of the present invention, the quasi- 
equilibrium nuclear spin ensemble state lasts for a time of at least 3T L Thus, the method 
may further comprises allowing said quasi equilibrium nuclear spin ensemble state to remain 
for a time of equal to or greater than 3Ti, where Ti is the spin lattice relaxation time before 
breaking at least one symmetry operation of said Hamiltonian; 

[001 9] In any of the above aspects of the present invention, a sample is prepared 
which has a Hamiltonian possessing one or more symmetry operations, this state can exist 
for a time of at least 3Ti before one or more of the symmetry operations is removed from the 
Hamiltonian. 

[0020] A number of different methods may be used to impose and remove a 
symmetry operation for the nuclear spin Hamiltonian operator and some of these are 
dependent on the nature of the sample. 

[002 1 ] When the molecules of the sample comprise nuclear spins of different 
isotopic types or have different chemical shifts and the molecules have sufficient geometrical 
symmetry said step of providing a sample where the nuclear spin Hamiltonian operator of the 
component molecules of the sample possess one or more symmetry operations may comprise 
reducing a magnetic field to a value where a symmetry operation is imposed on said 
Hamiltonian. A molecule of sufficient geometrical symmetry is defined as a molecule for 
which the spin Hamiltonian, including all direct and indirect spin-spin couplings, has higher 
symmetry when the chemical shifts are ignored than when the chemical shifts are present. 
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[0022] Alternatively, the sample could be kept in high field while a radiofrequency 
pulse sequence is applied which has the effect of suppressing the chemical shift differences, 
see for example J. D. Ellett and J. S. Waugh, J. Chem. Phys. 51, 2851 (1969); M. H. Levitt, 
R. Freeman and T. A. Frenkiel, Adv. Magn. Reson. 1 1, 47 (1983). 

[0023] Alternatively, a sample may be subjected to a chemical reaction, configured 
to switch the molecules between different conformational or configurational states, in which 
spin symmetry operations are imposed in order to provide a sample where the nuclear spin 
Hamiltonian operator of the component molecules of the sample possess one or more 
symmetry operations. An example of such a chemical reaction is a photoisomerization 
reaction. Such a reaction may be used to change a molecule from a mixed cis-trans form to a 
pure cis or pure trans form. 

[0024] A number of different methods may be employed for creating the quasi- 
equilibrium nuclear spin ensemble state. For example, the sample may be subjected to 
radiofrequency pulses of magnetic field in a high magnetic field environment. Specifically, 
thermal equilibration of a sample in high magnetic field, perturbation of the equilibrium of 
the nuclear spin system using applied radiofrequency pulses, and transport of the sample into 
the low magnetic field. 

[002 5] Those skilled in the art are aware of the technique of manipulating the 
populations of spin states using timed sequences of pulsed radiofrequency fields, resonant 
with the precession frequency of the nuclei in the high static magnetic field. The sequence 
of pulsed fields may be configured to produce a non-equilbrium state in high field, which 
gives rise to a quasi-equilibrium configuration of spin state populations, once the sample is 
transported to low magnetic field. 

[0026] The quasi-equilibrium nuclear spin ensemble state may also be provided by 
reacting an unsaturated symmetric molecule with parahydrogen. 

[002 7] Alternatively, the quasi-equilibrium nuclear spin ensemble state may also be 
provided by reacting an unsaturated non-symmetric molecule with parahydrogen, the product 
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of said reaction being placed in a magnetic field which is sufficiently low such that the 
difference in chemical shifts be ignored. 

[0028] Providing a sample with the desired symmetry and creating a quasi- 
equilibrium nuclear spin ensemble state, may be performed as distinct steps in sequence or 
otherwise, or may be performed together, for example by manipulating the populations of the 
spin states of a molecule such that a quasi-equilibrium spin state is formed once a symmetry 
operation is provided to the Hamiltonian. 

[0029] Since there are no relaxation processes internal to the molecule which can 
convert the quasi-equilibrium nuclear spin ensemble state to a fully equilibrated ensemble 
state, the lifetime of the quasi-equilibrium nuclear spin ensemble state can be much longer 
than the conventional relaxation time Ti.. For this reason, such states may be used a long- 
term repository of nuclear spin order. For example, the inventors found that, in systems of 
inequivalent spin pairs in low magnetic field, the decay of the singlet state population, or 
other quasi-equilibrium nuclear spin ensemble state, is governed by a time constant which is 
much longer than the spin relaxation time T\ This long lifetime allows nuclear spin order to 
persist during processes which take longer than Tj. 

[0030] Thus, while said quasi-equilibrium nuclear spin ensemble state is maintained, 
the sample may flow through an apparatus or a body, may be reacted with a second sample, 
may be purified or placed in contact with a second sample, in order to allow a transfer of 
nuclear spin order. The quasi-equilibrium nuclear spin ensemble state is preferably 
maintained for at least 5Ti, more preferably at least lOTj. 

[003 1 ] The quasi-equilibrium nuclear spin ensemble state may be observed by (i) 
breaking the symmetry operation of the molecular spin Hamiltonian followed by (ii) a 
sequence of events leading to the generation of an informative NMR signal. 
[0032] The operations leading to the breaking of symmetry are the same as those 
discussed above in the context of the making of symmetry. For example, in systems 
displaying chemical inequivalence, the symmetry may be broken by imposing a sufficiently 
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large external magnetic field. In other cases, the symmetry may be broken by 
photoisomerization processes or other chemical reactions. 

[003 3] Applying a sequence of magnetic fields to generate an NMR signal from said 
sample may be performed in either a high static magnetic field background or a low static 
magnetic field background. A low field is considered to be a field which is equal to or less 
than 20mT. 

[0034] The condition for a low field is that the difference in the precession 
frequencies of the involved spins must be equal to, or smaller than the J-couplings. In the 
case of two protons, with a chemical shift difference of lOppm and a J-coupling of 1 Hz 
(roughly worst case), this condition implies that the field should be less than around 2mT, 
which is easy to achieve, even without magnetic shielding. If the chemical shift difference is 
smaller, or if the J-coupling is larger, larger magnetic fields may be tolerated, in direct 
proportion to the J-coupling, and in inverse proportion to the chemical shift difference. The 
heteronuclear case is more difficult; for example, for a system of protons and carbon- 13 
nuclei, with a J-coupling of 120 Hz, the external field must be less than 3 pT, which requires 
careful magnetic shielding. A high field condition as used herein is when the difference in 
the precession frequencies of the involved spins is larger than the J-couplings. 
[003 5] The NMR signal may be detected in the standard way using electromagnetic 
induction within a coil. Alternatively, a direct magnetic field detecting device such as 
Superconducting Quantum Inteference Device "SQUID" may be used (see for example R. 
McDermott, A. H. Trabesinger, M. Muck, E. L. Hahn, A. Pines and J. Clarke, Science 295, 
2247-2249 (2002)). 

[0036] In a third aspect, the present invention provides an NMR apparatus, the 
apparatus comprising: means for providing a sample where the nuclear spin Hamiltonian 
operator of the component molecules of the sample possess one or more symmetry 
operations; means for creating a quasi equilibrium nuclear spin ensemble state in a sample, 
said quasi equilibrium nuclear spin ensemble state comprising at least two manifolds of spin 
states which transform differently under said symmetry operations of said Hamiltonian and 
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said manifolds having different mean nuclear spin populations, said quasi equilibrium 
nuclear spin ensemble state not having a singlet state population; means for breaking the 
symmetry operation of said Hamiltonian; means for applying a sequence of magnetic fields 
to generate a nuclear magnetic resonance signal from said sample; and means for detecting 
said nuclear magnetic resonance signal. 

[003 7] The means for providing a sample where the nuclear spin Hamiltonian 
operator of the component molecules of the sample possess one or more symmetry 
operations may be provided by a chemical or physical environment configured to produce 
the desired sample. For example, a means may be provided to perform any of the above 
suggested methods for providing such a sample. 

[0038] Similarly, means for creating a a quasi equilibrium nuclear spin ensemble 
state in said sample may be a means for allowing performance of any of the above described 
methods for creating the a quasi equilibrium nuclear spin ensemble state. 
[003 9] Although the above discussion has concentrated on an NMR method, the 
creation of a quasi equilibrium nuclear spin ensemble state may be used for other 
applications, for example for storing information in the form of quantum bits or "qubits". 
[0040] Thus, in a fourth aspect, the present invention provides a method of storing 
information, said method comprising: providing a sample where the nuclear spin 
Hamiltonian operator of the component molecules of the sample possess one or more 
symmetry operations; and creating a quasi equilibrium nuclear spin ensemble state in a 
sample, said quasi equilibrium nuclear spin ensemble state comprising at least two manifolds 
of spin states which transform differently under said symmetry operations of said 
Hamiltonian and said manifolds having different mean nuclear spin populations; assigning a 
bit value to at least one spin state of said quasi equilibrium nuclear spin ensemble state and 
storing said spin state for a time equal to or substantially greater than 3Ti, where Ti is the 
spin lattice relaxation time. 

[004 1 ] In a fifth aspect, the present invention provides a method of storing 
information, said method comprising: providing a sample where the nuclear spin 
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Hamiltonian operator of the component molecules of the sample possess one or more 
symmetry operations; and creating a quasi equilibrium nuclear spin ensemble state in a 
sample, said quasi equilibrium nuclear spin ensemble state comprising at least two manifolds 
of spin states which transform differently under said symmetry operations of said 
Hamiltonian and said manifolds having different mean nuclear spin populations, said quasi 
equilibrium nuclear spin ensemble state not having a singlet state population; assigning a bit 
value to at least one spin state of said quasi equilibrium nuclear spin ensemble state and 
storing said spin state for a time equal to or substantially greater than 3Ti, where Ti is the 
spin lattice relaxation time. 

[0042] The preferred methods discussed above in relation to NMR for providing the 
sample and creating the quasi equilibrium nuclear spin ensemble state may also be used with 
the above fourth and fifth aspects of the invention. 

[0043] The above methods may be used to store qubits in quantum computing or 
related applications. Previously, NMR techniques have been suggested for quantum 
computation by: D. G. Cory, A. F. Fahmy and T. F. Havel, Proc. Natl. Acad. Sci. USA 94, 
1634-1639 (1997); N. A. Gershenfeld and I. L. Chuang, Science 275, 350-356 (1997). The 
techniques covered in these papers may be combined with the present invention. 
[0044] In a sixth aspect, the present invention provides an apparatus for storing 
information, said apparatus comprising: means for providing a sample where the nuclear spin 
Hamiltonian operator of the component molecules of the sample possess one or more 
symmetry operations; and means for creating a quasi equilibrium nuclear spin ensemble state 
in a sample, said quasi equilibrium nuclear spin ensemble state comprising at least two 
manifolds of spin states which transform differently under said symmetry operations of said 
Hamiltonian and said manifolds having different mean nuclear spin populations; means for 
assigning a bit value to at least one spin state of said quasi equilibrium nuclear spin ensemble 
state and storing said spin state for a time equal to or substantially greater than 3Ti, where Tj 
is the spin lattice relaxation time. 
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[0045] In a seventh aspect, the present invention provides an apparatus for storing 
information, said apparatus comprising: means for providing a sample where the nuclear spin 
Hamiltonian operator of the component molecules of the sample possess one or more 
symmetry operations; and means for creating a quasi equilibrium nuclear spin ensemble state 
in a sample, said quasi equilibrium nuclear spin ensemble state comprising at least two 
manifolds of spin states which transform differently under said symmetry operations of said 
Hamiltonian and said manifolds having different mean nuclear spin populations, said quasi 
equilibrium nuclear spin ensemble state not having a singlet state population; means for 
assigning a bit value to at least one spin state of said quasi equilibrium nuclear spin ensemble 
state. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0046] The present invention will now be described with reference to the following 
preferred non-limiting embodiments in which: 

[0047] Figure 1 is a flow diagram outlining a method in accordance with an 
embodiment of the present invention; 

[0048] Figure 2 is a flow diagram to illustrate the switching of symmetry operations 
in the Hamiltonian; 

[0049] Figure 3 is a flow diagram illustrating a method in accordance with an 
embodiment of the present invention; 

[0050] Figure 4a is a schematic plot of B-field applied to the sample against time.; 
[005 1 ] Figure 4b is a schematic plot of radiofrequency magnetic pulse sequences 
applied perpendicular to the sample against time; 

[0052] Figure 4c is a schematic plot of the distribution of spin state populations of a 
sample subjected to the static fields of Figure 4a and the pulse sequences of Figure 4b; 
[005 3] Figure 5a is a conventional ! H NMR spectrum of the solution of 2,3- 
dibromothiophene (inset) in DMSO-d 6 ; 

[0054] Figure 5b is a spectrum generated by the method of Figures 3 and 4a; 
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[005 5] Figure 5c is a spectrum generated by the method of Figures 3 and 4a, using a 
different pulse sequence to t used to produce the spectrum of Figure 5b; 
[0056] Figure 6 schematically illustrates a plot of signal amplitude against time when 
the sample is placed in a low field condition; 

[0057] Figures 7a through 7e are plots of singlet population against time for a 
number of different magnetic fields which are to be used as the low field, wherein Figure 7a 
is for a low field of 0T, Figure 7b is for a low field 20mT, Figure 7c is for a low field of 
200mT, Figure 7d is for a low field of 500mT, Figure 7e is for a low field of IT and Figure 
7f is for a field of 5T; and 

[0058] Figure 8 is a flow diagram schematically illustrating a method in accordance 
with a further embodiment of the present invention; 

[0059] Figure 9 is a diagram schematically illustrating a chemical reaction of a 
symmetric molecule with parahydrogen; 

[0060] Figure 10 is a flow diagram schematically illustrating a method in accordance 
with a further embodiment of the present invention; 

[006 1 ] Figure 1 1 is a diagram schematically illustrating a chemical reaction of a non- 
symmetric molecule with parahydrogen; 

[0062] Figure 12 is a flow diagram schematically illustrating a method in accordance 
with a further embodiment of the present invention; 

[0063] Figures 13a and 13b schematically illustrate a photoisomerization process 
which may be used to impose or remove a symmetry operation on a Hamiltonian; and 
[0064] Figure 14 is a flow diagram schematically illustrating a method in accordance 
with a further embodiment of the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 
[0065] Figure 1 outlines the basic steps in accordance with an embodiment of the 
present invention. 
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[0066] In step S101, a sample is provided where the nuclear spin Hamiltonian 
operator of the component molecules of the sample possess one or more symmetry 
operations. This step may involve providing a sample which has this quality inherently 
(such a method will be described in more detail with reference to figure 8), alternatively, the 
sample may only possess this quality in certain environments (this will be described with 
reference to figures 3, 10 and 14), as a further option, the sample may need to be chemically 
modified to achieve the symmetry requirement (this will be described with reference to 
figure 13). 

[0067] In Step SI 03, a quasi-equilibrium nuclear spin ensemble state is created, how 
this state may be created will be described later with reference to figures 3 to 14. The term 
quasi equilibrium nuclear spin ensemble state is used to refer to the state of a spin ensemble 
in which manifolds of different symmetry have reached internal equilibrium but in which the 
different manifolds are not equilibrated with each other. 

[0068] Once the quasi equilibrium nuclear spin ensemble state has been created, it 
may be stored for a reasonably long time, much longer the spin-lattice relaxation time Tj. Tj 
is believed to be the time constant which determines when the spin memory is quenched in 
NMR experiments. However, the inventors have found that the quenching of quasi 
equilibrium nuclear spin ensemble states is governed by a time constant T qe which is much 
greater than Tj. This allows the sample to undergo processes SI 05 before the quasi 
equilibrium nuclear spin ensemble state decays and for a time greater than or equal to 3Ti. 
For example, the sample may be injected or inserted into a human or animal body and its 
passage through the body may then be studied using NMR. Since the sample may be 
allowed to propagate in the body for longer than conventionally allowed due to the discovery 
of the long lifetime of the quasi equilibrium nuclear spin ensemble state, more information 
may be gathered and different types of experiments may be conducted utilising the long 
lifetime of the quasi equilibrium nuclear spin ensemble state. 

[0069] The extension of spin memory lifetimes beyond Tj by using quasi-equilibrium 
states is distinct from the extension of spin dephasing lifetimes beyond T 2 demonstrated in 
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the context of Decoherence-Free Subspaces as employed in quantum computation (L. Viola, 
E. M. Fortunato, M. A. Pravia, E. Knill, R. Laflamme and D. G. Cory, Science 293, 2059- 
2063 (2001); E. M. Fortunato, L. Viola, M. A. Pravia, E. Knill, R. Laflamme, T. F. Havel 
and D. G. Cory, quant-ph /0210057 (2002); D. A. Lidar and K. B. Whaley, quant-ph 
/0301032 (2003); J. E. Ollerenshaw, D. A. Lidar and L. E. Kay, quant-ph /0302175 (2003)). 
The main distinctions are: (i) the major Ti mechanism is usually the modulation of the 
dipole-dipole coupling by random molecular motion, while the experiments reported in the 
quantum computation field demonstrate the compensation of an artificial shortening of T2 by 
molecular diffusion in strong magnetic field gradients; (ii) the methods reported in the above 
references are incapable of demonstrating a lengthening of the spin memory time beyond Ti. 
[0070] The sample may be reacted with other samples or may undergo a purification 
process during step SI 05. Alternatively, the sample may be placed in physical contact with 
other molecules to transfer polarisation of the quasi equilibrium nuclear spin ensemble state. 
[0071 ] The quasi equilibrium nuclear spin ensemble state may only be observed if 
the symmetry of the sample is broken SI 07. This may be done by a number of methods 
which will be described with reference to figures 2 to 13. 

[0072] The sample is then subjected to a sequence of magnetic fields to generate 
nuclear spin magnetization. The sequence of magnetic fields may be applied in a high 
magnetic field environment as described in detail with reference to Figure 4, or a low or zero 
magnetic field environment as described in R. McDermott, A. H. Trabesinger, M. Muck, E. 
L. Hahn, A. Pines and J. Clarke, Science 295, 2247-2249 (2002). 

[0073] The signal from the NMR active sample can then be detected using either a 
high field detection technique such as electromagnetic induction of current in a coil or direct 
detection of the magnetisation using a SQUID or the like. 

[0074] Figure 2 schematically illustrates the symmetry considerations for steps SI 03 
to SI 09 of figure 1. The small rectangles represent two different state manifolds, and the 
shading represents their populations. In S201, the quasi equilibrium nuclear spin ensemble 
state has been prepared and the sample has a Hamiltonian which has one or more symmetry 
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operations. In this environment, the quasi equilibrium nuclear spin ensemble state decays 
with a long time constant T qe , which is much longer than the spin lattice relaxation time T|. 
[0075] In step 203, the symmetry operations are removed from the Hamiltonian, 
allowing the NMR signal to be observed,. Due to spin-lattice relaxation, the spin population 
between the states equilibrates and an equilibrium state is reached after Tj in step S205. 
[0076] Figure 3 and figure 4 schematically illustrate a method in accordance with an 
embodiment of the present invention. Figure 3 is a flow diagram summarising the processes, 
while Figures 4a, 4b and 4c describe the changes to the sample and its environment. 
[0077] In NMR procedures, a sample is subjected to a high constant magnetic field 
Bhigh. A much smaller rf pulsed magnetic field is then provided perpendicular to the high 
magnetic field Bhigh- 

[0078] The pulsed magnetic field follows a pulse sequence A+. The pulse sequence 
is applied at the Larmor frequency of the spins in the high magnetic field. The pulse 
sequence will be described in greater detail later in the description. It is well known to those 
skilled in the art of NMR how to design pulse sequences in order to achieve a desired 
distribution of spin states within a sample. For the purposes of symmetry-breaking, the high 
magnetic field simply needs to be an order of magnitude larger than the "low" field discussed 
above, i.e. more than a few mT for typical proton systems, and more than a few nT for 
heteronuclear systems. However, for the purposes of sensitive observation of the NMR 
signals, the field will usually be as large as possible, i.e. between 1 and 20 Tesla, depending 
upon the instrument available. 

[0079] Figure 4c schematically illustrates how the spin states may be manipulated 
using pulse sequence A+. 

[0080] At time point 1, the nuclear spin system has four Zeeman energy levels with 
energies in the order E aa (3) < E a p (5) < Ep a (7) < Epp (9). In the figure, E a p (5) and Ep a (7) 
are shown to be virtually identical. However, there is a slight energy difference between 
these two levels. The above assumes that the chemical shift 62 is greater than the chemical 
shift 81, the gyromagnetic ratio yand the J-coupling are both positive and that the spin 
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system is weakly-coupled in high fields. The population is distributed throughout the four 
states. In thermal equilibrium at ambient temperature, the lower most level | aa> has a 
small excess population (depicted by the filled balls), while the uppermost level | PP> has a 
slightly depleted population (depicted by the white balls). These population differences are 
very small in the case of conventional thermal polarisation, but may be of the order of one in 
spin hyperpolarization experiments. The sample is in thermal equilibrium and this stage of 
the sample corresponds to step S301 in figure 4. 

[0081 ] The pulse sequence A+ is then applied until time point 2, taking the sample to 
step S303 in figure 4. During this time, the population of state | otoc> exchanges with that of 
| otp>, and that of | (3(X> with that of | pp>. Thus, a non-equilibrium population distribution 
is generated, in which level | ocp> has an excess population, while | poo has a depleted 
population. This is shown at time point 2. 

[0082] During the interval ^ from time point 2 to time point 3, the sample is 
transported to a region of low magnetic field Bi ow , which is of the order of milliTesla, taking 
the sample to S305 of figure 4. Typically, Bi ow will be at most 20mT. In low field, it is 
convenient to discuss the behaviour of the nuclear spins using the four eigenstates of the pure 
J-coupling Hamiltonian [1]. 

Hj=2nJI r I 2 [1] 

[0083] The four J-coupling eigenstates may be classified as the three components of 
a triplet state, plus a singlet state. The triplet state components are denoted by equation [2]. 

\T { )=\aa) 

|r 0 )=-LM+lA*) PI 

[0084] The triplet states have the same J-coupling eigenvalue of + \ 7tJ '. The singlet state 
is denoted by equation [3]. 
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\S 0 )=j=$af})-\/Ja)) [3] 

[0085] The singlet state has a J-coupling eigenvalue of -\nJ . The triplet and singlet 

states are energy eigenstates of the nuclear spins only if the magnetic field is zero. The 
energy level diagram in this case is depicted next to time point 3 in Figure 4c. Although we 
use the zero-field eigenstates as a basis for a representation of the spin density operator, the 
treatment below does not assume a field of exactly zero. 

[0086] If the transport from high to low field is fast compared to T\ but slow compared 
to the J-coupling, the excess population of the high-field eigenstate | ap> is transferred 
adiabatically into an excess population of the low-field eigenstate | So>. Similarly, the 
depleted population of the high-field eigenstate | (3a> is transferred adiabatically into a 
depleted population of the zero-field eigenstate | To>. The idealized distribution of the 
singlet and triplet state populations is shown at time point 3 in Figure 4c. This assumes pure 
adiabatic transfer and neglects Tj relaxation during the transport interval. 
[0087] The nuclear spins evolve in low magnetic field under a combination of coherent 
and incoherent relaxation processes. Here we use a standard second-order perturbation 
treatment of the modulate dipole-dipole interaction to analyze the spin-lattice relaxation in 
low field. 

[0088] In general, the spin density operator p evolves according to the Liouville-von 
Neumann equation [4]. 

at 

[0089] Where the low-field Liouvillian is given by equation [5]. 

L L F = ^LF~ MlF [5] 

[0090] H LF is the superoperator for the commutator with the low-field spin Hamiltonian 
of equation [6]. 

H LF = -yB LOW (1 + S x )/ lz - }B LOW (1 + S 2 )I 2Z + Hj [6] 
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[009 1 ] The thermal polarization terms may be ignored in low field. If the relaxation is 
dominated by the intramolecular dipole-dipole coupling mechanism, the relaxation 
superoperator may be evaluated by second-order time-dependent perturbation theory and is 
given by equation [7]. 

t LF =-^A DD {-T)A DD (0)dT [7] 

[0092] H DD (t) is the commutation superoperator for the dipole-dipole coupling 
Hamiltonian, evaluated at time t, and the overbar represents an ensemble average. Equation 
[4] is general for any field strength low enough to satisfy the extreme narrowing condition 
for the dipole-dipole relaxation. 

[0093] It is convenient to us the matrix representations of the Liouvillian superoperator 

A 
A 

L LF in a space of orthonormal spin operators. The following set of six operators generate a 
suitable basis: 



11) 


ISoXSol 


= ±(i?i? + i?i a 2 -i;r 2 -i;r 2 ) 


12) 


|T,XT,| 




13) 


IToXTol 




14) 


|T-i><T.,| 


=lflf 


15) 


=^S 0 )(T 0 \ + \T 0 )(S 0 \) 


=-L(l?I?-lfl°) 


16) 


=iPo)(To\-\T 0 XS 0 \) 





[0094] where Ij are single-spin shift operators and \IJ 9 lf /are single-spin polarization 

operators. The ket 1 1) describes the population operators of the singlet state. Kets |2), |3) 
and 1 4) describe the population operators of the three triplet states. Kets | 5) and 1 6) 
represent the x and y-components of the coherence between the singlet state and the central 
component of the triplet. All of these operators commute with the angular momentum 
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operator along the field axis, as does the low-field Hamiltonian H LF . The operator basis of 

equation [8] is closed under the application of the low-field Liouvillian L LF . 

[0095] The matrix representation of the low-field relaxation superoperator T LF may be 
evaluated by standard relaxation theory, assuming pure intramolecular dipole-dipole 
relaxation driven by rigid isotropic random rotation with a correlation time x c . The low-field 
Liouvillian evaluates to 





f o 


0 


0 


0 


0 


-2 y >Q) l ° w 




0 


10 u L c 




\b\ 


0 


0 




0 




~\b 2 r c 


^b\ 


0 


2 Yl (O l ™ 




0 


\b\ 




10° l c 


0 


0 




0 


0 


0 


0 


~\b\ 


2nJ 






0 




0 


-2nJ 


-\b 2 r c 



[8] 

[0096] where the dipole-dipole coupling between the spins is given by b= - 
{^Anyfhf 3 , and r is the internuclear distance. The low-field chemical shift frequency 
difference is defined by equation [9]. 

< = -lS Um {d 2 -5) [9] 
[0097] In zero field co l ™ = 0 , the first row and column of equation [8] contains only 
zeros. This indicates that the singlet state population is dynamically isolated in low field, 
and is conserved under the evolution process. The zero-field singlet state may be used as a 
long-term repository for nuclear spin order, protected from intramolecular dipole-dipole 
relaxation. 

[0098] If the magnetic field is not exactly zero, equation [8] indicates that the singlet 
state is relaxed indirectly by dynamic couplings to the central triplet state and to singlet- 
triplet coherence. 

[0099] The idealized behaviour of the spin ensemble after an interval t L f in low field 
is depicted by the diagram at time point 4 in Figure 4c. If the field is sufficiently low, and 
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Tlf is long enough, the triplet populations equilibrate, while the singlet population is 
conserved. Thus, a quasi-equilibrium ensemble spin state is rapidly established. This quasi- 
equilibrium ensemble spin state is stable for a time much longer than T|. 
[001 00] The singlet population may be read out be adiabatic transport of the sample 
back into high field during the interval tf] nsp . This process breaks the symmetry of the 
sample and forms a sample which is an active NMR sample as illustrated at step S307 of 
figure 3. 

[001 01 ] Specifically, this process transforms the low-field singlet population into a 
population of the high-field | ocp> state, as depicted for time point 5 in Figure 4c. The 
population of the | ap> state is converted into observable NMR signals by the pulse 
sequence B at step S309 of figure 3. 

[001 02] Sequence B is given by three rf pulses separated by two unequal delays, and 
converts the | a(3> state population into antiphase NMR signals centered at the chemical 
shift 82. NMR signals deriving from Ti relaxation during the transport interval, on the other 
hand, only generate NMR signals at the chemical shift 81. The interesting NMR signals 
derived from low-field singlet storage are therefore cleanly separated from signals with a 
trivial origin. 

[001 03] The pulse sequence A+ is given by the two-pulse sequence 90o - Xi - 90_9o, 
where the symbol denotes a strong, non-selective pulse with flip angle |3 and phase § 
(both angles are specified in degrees, and the phases taken into account the sign of the 
precession and the radio-frequency mixing scheme. The delay Xi is set to the value 



For the experimental results, the delay was set to the value Xi = 2.05ms. 

[001 04] Standard spin operator theory shows that in the case 82 > 5i and y > 0, the pulse 

sequence A+ has the approximate propagator 




where 



[10] 
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U(A + )= exp[-/|(/ lz +/ 2Z )|exp{+/^} [11] 



neglecting the effect of the J-coupling during the interval Xi. This assumption is valid in the 
case that the chemical shift frequency difference is much larger than the J-coupling (weak 
coupling approximation). The transformations of the individual state kets are 

U(A + ) | aoo = - 1 ap> 

U(A + ) | poo = -i | aa> 

£/(A + )|ap> = -i|pp> 

C/(A + )|pp>- | poo 

which leads to the exchanges of populations sketched in Figure 3c. The population of the 
I aP> state at time point 2 is derived from the thermal equilibrium population of the | aoo 
state. 

[00 105] The pulse sequence B is given by 90o - %i - 1 8O90 - Xfc - 9045, with pulse 



sequence delays given by T 2 = 



fy»A gh ) + X 4j ) and ^ = \y^ j \ ' F ° r the experimentai 



results, the delays were set to the values X2 = 44.64ms and X3 = 43.61ms. This sequence is 
designed to separate cleanly the interesting NMR signals deriving from the low- field singlet 
state from the uninteresting signals generated by relaxation of the spin system during the 
second transport interval. 

[00 1 06] From standard spin-operator theory, the propagator U(B) for this pulse sequence 
transforms the population operator of state | ocp> into a superposition of many operator 
terms, including some that represent antiphase single-quantum coherences: 

[001 07] The (-l)-quantum coherences represented by operators and ifl^ generate 
the diagnostic antiphase signals near the chemical shift 02. 



4 [12] 
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[00 1 08] Sequence B gives rise to no net signals at the chemical shift 82 when applied to 
in-phase Zeeman magnetization, created by Ti relaxation during the second transport interval 
T ?2s P - This may be seen from the transformations 

U(By iz U(Bf = -± expf / */ A p« - 1 expf/ % };i° 2 
1 1 [13] 

+\ ex pf ■ % y*t + \ ex pf ■ % y*i 

and 

u(B)i 2Z u(Bf = -\n+i\n-i\i;r 2 +|/,-/ 2 - m 

[00 1 09] Partially-relaxed Zeeman magnetization of spin Ii only generates signals at the 
chemical shift 81 , while partially-relaxed Zeeman magnetization of spin I 2 only generates 
multiple-quantum coherences, which do not generate a signal at all. In principle, cross- 
correlated relaxation during the transport interval could also generate spurious signals at the 
shift 82. 

[001 1 0] Figure 5 illustrates NMR spectra for a sample of 500 microlitres of a 17.7 mM 
solution of 2,3-dibromothiophene in DMSO-d 6 at 20°C. The sample was contained in a 
5 mm high-resolution NMR tube equipped with a Young valve to facilitate degassing. The 
sample contained a water impurity in approximately 27 mM concentration. The sample was 
subjected to three freeze-pump-thaw degassing cycles, each lasting over 30 minutes, to 
remove dissolved oxygen. 

[001 1 1 ] The NMR spectra illustrated in Figure 5 were obtained using a Varian Infinity + 
400 MHz NMR system at a field of B h i g h * 9.4 T, using an actively shielded 89 mm bore 
magnet. A standard 5 mm high resolution NMR probe was used. The 90° pulse duration 
was around 5 jlls. 

[001 1 2] The sample was transported between the high and the low fields using the 
standard pneumatic sample elevator. The pressure of the elevator was increased during the 
lift operation, and reduced during the insert operation in order to accelerate the transport 
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times. Reproducible transport times of r,;^ = 9±lsand = 13 ± 2s measured from 

when the sample is fully loaded into the high field generating coil to when the sample is 
outside the magnet bore. 

[001 1 3] The field Bi ow was estimate using a Hall-effect Gauss meter, with zero field 
calibrated to an accuracy of approximately 10 using a mu-metal shield far from the 
magnet. 

[001 1 4] For a low or storage field of Bi ow = 20 mT, the sample was elevated to the top of 
the magnet bore and allowed to remain there for the time Tlf- For even lower storage fields, 
the sample was carried physically to calibrated locations within the stray field of the magnet, 
or in the case Bi ow = 0, a magnetically shielded chamber. 

[001 1 5] In the above sample, the proton spin system is defined by 8i= 7.1 1 ppm, 82= 7.72 
ppm and J = 5.7 Hz. Figure 5a illustrates a conventional proton NMR spectrum, obtained at 
400 MHz in the field of B h igh « 9.4 T. The two doublets of a typical AX spin system 51, 53 
are clearly visible. 

[001 1 6] A conventional inversion recovery series (not shown), demonstrates that the spin- 
lattice relaxation constants for the inequivalent sites are very similar in high field (Ti(5i) = 
16.5±0.2s;T,(62)=17.1±0.2s). 

[001 1 7] Figure 5b is a plot of the NMR spectrum for the above sample after it has been 
subjected to the method described with reference to figures 3 and 4. A quasi-equilibrium 
state was formed by subjected the sample to a high magnetic field and pulse sequence A+, 
then transferring the sample adiabatically back to a low storage field. The spectrum in figure 
3b was taken after the sample was maintained in a low field for a time Xlf of 100s (i.e. much 
longer than Ti for either 8i or 82). The spectrum of figure 5b shows diagnostic antiphase 
signals 55 at around 52=7.72 ppm, and an antiphase dispersion structure 57 around 81 = 7.1 1 
ppm. Both of these structures are due to the creation of the quasi-equilibrium state in the 
preparation of the sample. 

[001 1 8] The inset of figure 5b shows the region from 7.76 to 7.68 expanded. 
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[00 1 1 9] The quasi-equilibrium state in the sample is created by manipulating the spins 
using pulse sequence A+ as previously described. To prove that the creation of the singlet 
state is responsible for the features which exist long after Ti, a second experiment was 
performed. 

[001 20] The second experiment was performed in the same manner as the first 
experiment, but used the pulse sequence A_ instead of the pulse sequence A+. Pulse 
sequence A. is derived from A+ by changing one of the pulse phases. This exchanges the 
population of state |ap> with that of |pp> instead of with |otot>. The population of |ap> is 
therefore depleted at time point 2 as opposed to enhanced. 

[001 21] In detail, pulse sequence A. is given by two pulse sequences 90i 8 o-Xi-90. 90 . The 
sequence differs from A+ by a 1 80° phase shift of the first pulse. The spin propagator for 
this sequence, under the weak coupling approximation, is: 

U(A_)= e xp[-/|(/ 12+ / 2Z )|exp{ + /^} [15] 

which leads to the following ket transformations: 

U(A.) | aoo = - 1 ap> 
t/(A.)lpa> = i|pp> 
<7(A_) | ap> = i | aoo 
C/(A.)IPP>= I pa> 

[001 22] The population of the | aP> state at time point 2 is derived from the thermal 
equilibrium population of the | aoo state. 

[001 23] Since the | aP> state is connected adiabatically to the low field singlet, the 
sequence A. as well as A+ can be used to prepare a low field singlet. However, the 
population of the low field singlet deviates from the mean in opposite senses for the states 
formed using A+ and A.. This variation from the mean in opposite senses means that the 
diagnostic antiphase signals 55 will change sign. 

[001 24] Thus, in figure 5c, the diagnostic antiphase signals 59 at 62 = 7.72 ppm have a 
change in sign from the equivalent signals 57 in figure 2b. The signals around 8i=7.1 1 ppm 
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are also perturbed, but in more complicated way since these signals are superpositions of 
several contributions. 

[001 25] Since the total interval ^L^+^f+^U = 1205 exceeds Ti by a factor of seven, 
the spectra of figures 5b and 5c prove that the memory of the nuclear spin system has been 
extended significantly beyond the T| limit. 

[001 26] Figure 6 is a plot of the amplitude of the 7.72 ppm antiphase signals against as 
function of x L f for storage magnetic field strengths of 0T, 2mT and 20 mT. The amplitude of 
the 7.72 ppm antiphase signal is related to the singlet population. 
[00 127] In a storage field of Bi ow -l .8 ± 0.5 mT, the time constant for the quasi- 
equilibrium decay is estimated to be T qe = 104 ± 5s. The quasi-equilibrium decay time 
constant had no appreciable field dependence for magnetic fields less than Bi ow = 20 mT. 
The effect of higher magnetic fields on the population of the singlet state will be described 
with reference to figure 5. 

[001 28] The quasi-equilibrium decay is observed to accelerate when the proton 
concentration in the solution is increased. This indicates that intermolecular proton-proton 
relaxation is a strong contributor to the quasi-equilibrium relaxation, and that highly dilute 
solutions are advantageous to achieve long lifetimes. Quasi-equilibrium relaxation may also 
be caused by paramagnetic impurities such as dissolved oxygen, intermolecular interactions 
with solvent deuterons, and scalar relaxation of the second kind via the Br nuclei. 
[001 29] Figures 7a to 7f are plots of simulations of the singlet state population as a 
function of time for a plurality of different fields. In figure 7a Bi ow = 0T, figure 7b 
B] O w=20mT, figure 7c Bi ow =200mT, figure 5d Bi ow =500mT, figure 5e Bi 0W =lT and figure 5f 
Biow=5T. 

[00 1 30] The trajectory of the singlet population in a magnetic field is predicted using the 
Liouville-von Neumann equation [4]. For 2,3-dibromothiophene, the distance between 
proton nuclei was estimated to be r = 253 pm. The spin-lattice relaxation in high field was 
assumed to be dominated by intramolecular dipole-dipole relaxation, modelled as the 
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isotropic rotational diffusion of a rigid molecule. With this model, the observed Ti value of 
20s corresponds to a rotational correlation time of T c =15 ps. 

[00 131] At low values of Bi ow , the singlet state population is predicted to be time 
independent over a timescale of hundreds of seconds as seen in figures 7a and 7b. As the 
storage field is increased beyond around 200 mT, a rapid low amplitude oscillation of the 
singlet state population is imposed on a slower decay as can be seen in figures 7c and 7d. 
Finally, a higher fields, the amplitude of the oscillations becomes very large and at long 
times the population equilibrates among the four available states, indicating loss of all spin 
order as can be seen in figures 7e and 7f 

[001 32] The experimental results show a significant damping of the quasi-equilibrium 
spin order on the timescale of 100s for fields much less than 200mT. This indicates the 
participation of relaxation mechanisms other than intramolecular dipole-dipole coupling. 
Intermolecular dipole-dipole relaxation involving water protons or protons on different DBT 
molecules is likely to be the major quasi-equilibrium damping mechanism. 
[00 133] In the method described with reference to figures 3 and 4, symmetry operations 
are imposed on or removed from the Hamiltonian by moving the sample from a high to a low 
magnetic field. However, other options are possible. 

[00 1 34] Figure 8 is a flow diagram for performing a method in accordance with an 
embodiment of the present invention where the sample 401 is inherently symmetric in both a 
high or a low magnetic field. An example of such a sample 401 is shown in figure 9. 
[00135] A quasi equilibrium nuclear spin ensemble state may be formed by reacting a 
sample with parahydrogen in a high field environment, the "Pasadena" effect (Russell 
Bowers et al. J. Am Chem. Soc, 109 5541-5542 (1987)), or in a low field environment, the 
"Altadena" effect (Pravica et al, Chemical Physics Letters, 145, 255-258 (1988)). 
[001 36] The symmetric sample 401 is reacted with parahydrogen 403 the presence of a 
catalyst at step S405. This reaction is shown schematically in figure 9. Parahydrogen is a 
form of hydrogen gas in which the hydrogen nuclei are in the nuclear spin singlet state. In 
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figure 9, R refers to identical chemical groups. A sample having a quasi equilibrium nuclear 
spin ensemble state is produced 407. 

[001 37] The sample with the quasi equilibrium nuclear spin ensemble state 407 may then 
be subjected to processes for a time exceeding 3Ti as explained in relation to figure 1. 
[00 138] To break the symmetry of the sample and thus to allow an NMR signal to be 
observed, the sample is subjected to a photoisomerization process in step S409 by briefly 
illuminating the sample with light flash 411. 

[001 39] A non-equilibrium state of a non-symmetric product 413 is produced as a result 
of the photoisomerization process and the sample is investigated using NMR techniques in 
stepS415. 

[001 40] In the above described method, the sample is kept in a high field environment all 
the time and symmetry is switched on an off using either a chemical reaction or 
photoisomerization. 

[00 141] Figure 10 schematically illustrates a variation on the method described with 
reference to figure 8. Here, the starting sample 501 is an unsaturated non-symmetric 
molecule 501. An example of molecule 501 is shown in figure 11. 

[001 42] Molecule 501 is reacted with parahydrogen 503 in step S505 in the presence of a 
catalyst and in a low field environment to produce a sample 507 having a quasi equilibrium 
nuclear spin ensemble state. A quasi equilibrium nuclear spin ensemble state may be formed 
by introducing a sample to parahydrogen in a low field environment, the "Altadena" effect 
(Pravica et al, Chemical Physics Letters, 145, 255 to 258 (1988)). 

[001 43] In figure 1 1, A and B refer to non-identical chemical groups. In this case, a quasi 
equilibrium nuclear spin ensemble state is only produced if the external magnetic field is low 
enough that the chemical shift difference between the two hydrogen nuclei is negligible. 
[001 44] The symmetry is then broken by transferring the sample to a high field to produce 
non-equilibrium state 509. An NMR signal may then be detected in step S5 1 1 . 
[001 45] The flow chart of figure 12 schematically illustrates a variation of the method of 
figure 10. In the same manner as figure 10, unsaturated non-symmetric molecule 601 is 
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reacted with parahydrogen 603 in step S605 in the presence of a catalyst and in a high field 
environment to produce a sample 607 which is a non-symmetric hydrogenation product with 
enhanced spin order. 

[00 1 46] Next, a symmetry operation is imposed on the Hamiltonian of product 607 by 
transferring the sample adiabatically to a low field environment to form the quasi 
equilibrium nuclear spin ensemble state 609. 

[00 1 47] Once the sample 609 has been subjected to the desired processes, symmetry is 
broken by transferring the sample to a high field to produce NMR active sample 61 1 which 
can be measured in step S613. 

[00 148] As previously mentioned, symmetry operations may be imposed on or removed 
from the Hamiltonian using photoisomerization. 

[001 49] Chemical reactions, such as photoisomerization, may switch the molecules 
between different conformational or configurational states, in which spin symmetry 
operations are either imposed or removed. As an example, consider the molecule shown in 
figure 13 a. The two central CH2 groups are chemically equivalent, since the molecule has a 
mirror plane of symmetry. However, the double bonds are readily photoisomerized, giving 
the mixed cis-trans form shown in figure 13b. 

[001 50] The two central CH2 groups are now chemically inequivalent, since one has a cis 

ligand, while the other has a trans ligand. A further photoisomerization step generates the 

dX\-cis form, in which the two central CH 2 groups are chemically equivalent again: 

[001 5 1 ] Thus, a light flash 701 may therefore be used to switch the symmetry of the spin 

Hamiltonian. 

[001 52] The flow chart of figure 14 illustrates a further methodin accordance with an 
embodiment of the present invention. As in the case of themethod described with reference 
to figure 3, the sample is not symmetric when in thermal equilibrium in a high field 
environment 801. 

[001 53] A first RF pulse sequence is applied to the sample in a high field in order to 
create a non-equilibrium state 803. 
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[001 54] As before, the symmetry of the sample is removed, and due to the populations of 
non-equilibrium state 803, a quasi equilibrium nuclear spin ensemble state 805 is formed, 
where the sample 401 is inherently symmetric in both a high or a low magnetic field. An 
example of such a sample 401 is shown in figure 9. 

[001 55] In this case, the chemical shift differences between the nuclear spins are removed 
and symmetry is achieved by applying a special radiofrequency pulse sequence, a second rf 
pulse sequence, rather than transporting the sample into a low magnetic field. 
Radiofrequency pulse sequences with the appropriate properties are commonly used for 
heteronuclear spin decoupling in solution NMR (see for example, M. H. Levitt, R. Freeman 
and T. A. Frenkiel, Adv. Magn. Reson. 1 1, 47 (1983)). 

[001 56] The sample with the quasi equilibrium nuclear spin ensemble state 805 may then 
be subjected to processes for a time exceeding 3Ti as explained in relation to figure 1. 
Throughout this time, symmetry is maintained by the second rf pulse sequence which is 
continually applied to the sample. The sample is kept in a high field environment. 
[00 157] To break the symmetry of the sample and thus to allow an NMR signal to be 
observed, the second rf pulse sequence is removed. 

[001 58] A non-equilibrium state of a non-symmetric product 809 is produced as a 
result of removing the second rf pulse sequence. The sample is investigated using NMR 
techniques in step S81 1. 
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